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Switches in real systems take many forms, such as impacts, electronic relays, mitosis, and the 
implementation of decisions or control strategies. To understand what is lost, and what can be 
retained, when we model a switch as an instantaneous event, requires a consideration of so-called 
hidden terms. These are asymptotically vanishing outside the switch, but can be encoded in the form 
of nonlinear switching terms. A general expression for the switch can be developed in the form of a 
series of sigmoid functions. We review the key steps in extending the Filippov’s method of sliding 
modes to such systems. We show how even slight nonlinear effects can hugely alter the behaviour 
of an electronic control circuit, and lead to ‘hidden’ attractors inside the switching surface. 


There are increasingly general and powerful methods 
to study the dynamical effects of switches in otherwise 
smoothly evolving dynamical systems, whether a phase 
transition, a trade or policy decision, an electrical relay, 
a mechanical actuator or a biochemical valve. This paper 
is an attempt to highlight and deal with certain ambigu¬ 
ities of non-uniqueness of piecewise smooth systems like 
these. Attempts are often made to show the robustness 
and generality of discontinuous models by smoothing or 
regularizing a discontinuity, a process that is somewhat 
tautological because, as we show, entirely different dy¬ 
namics can be obtained with different choices of regular¬ 
ization. Fortunately this can be resolved using nonlinear 
switching terms, and we introduce a general method by 
developing sigmoid series expressions for discontinuous 
vector fields. We apply these ideas to a model of an el¬ 
ementary electronic control circuit, and give an example 
of how novel attractors may arise inside the switching 
surface. 

I. INTRODUCTION 

Many phenomena, in either our passive attempts to 
describe nature, or our active attempts to control it, 
are mixtures of steady behaviours and sharp transitions. 
The steady regimes are usually relatively easy to model. 
The transitions are often complicated and challenging to 
model, but occupy fleeting instants of time which make 
it tempting to simplify them, in the extreme, as discrete 
events in an otherwise smooth system. The effect of this, 
and how far we can develop the mathematics of a piece- 
wise determined system, is explored here. 

Modeling a switch by means of piecewise smooth func¬ 
tions is simple and convenient, requiring no detailed in- 
depth knowledge of the physical laws during switching, 
but it risks destroying information about the transition 
process. The resulting simplicity has inspired an exten¬ 
sive theory of discontinuity-induced bifurcations and sin¬ 
gularities. Much of the theoretical development is im¬ 
pressively powerful, and even highly rigorous (for reviews 
see ii iso, ill). Its generality, its applicability to real 


world switching phenomena, however, remains poorly un¬ 
derstood (see e.g. dll). We can, for example, derive 
piecewise smooth models as the limiting cases of smooth 
multiple timescale models, collapsing the fast switching 
dynamics into a single instant 12311 . or as describiM the 
average dynamics of stochastic |25j| or hysteretic [a, 
processes. But the limits are typically singular, and 
qualitatively different systems have the same piecewise 
smooth limit, resulting in ambiguity about the systems 
they describe. 

A resolution to this has emerged recently. It seems 
that generality in piecewise smooth systems requires a 
notion of hidden dynamics that is observable only during 
a switch and not outside it, encoded in nonlinear terms 
that represent the remnants of asymptotically small ef¬ 
fects that vanish in the limit of taking a piecewise smooth 
model. We shall present a few examples where nonlin¬ 
ear terms entirely alter the local and global behaviour of 
a system, including an application to electrical control 
systems. 

In engineering and biological applications particularly, 
an increasing number of examples of ‘nonlinear’ switches 
are becoming known. These are systems whose dynam¬ 
ical equation contain nonlinear combinations of steep or 
discontinuous ‘sigmoid’ functions, each sigmoid being an 
empirical representation of a particular switching action. 
In Boolean switching models of protein dynamics regu¬ 
lated via Hill functions, they are very common 
In [T^, it was shown that nonlinear switching can be 
used to model static friction without needing to augment 
Coulomb friction with a velocity dependence. In general, 
nonlinear switching can be used to model effects such as 
overshoot or lag, and can resolve deterministic ambigui¬ 
ties that arise when multiple switches coincide. It seems 
reasonable to hope that a deeper understanding of them 
will permit greater use for modeling across the physical, 
social, and biological sciences. 

In section m we recount the standard method for han¬ 
dling a discontinuity in a dynamical system, and give ex¬ 
amples of its ambiguity under smoothing. In section IIIII 
we discuss ambiguities in smoothing more generally. In 
section m we make these ambiguities explicit, and find 
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how to express them in the piecewise smooth system. We 
apply these ideas to a circuit model in section [V] In Sec¬ 
tion IVII we outline the method for solving a dynamical 
system at a discontinuity, followed by an example of the 
potential for novel dynamics, in the form of a hidden at¬ 
tractor, in section IVIIl A few closing remarks are made 
in section IVIIII 


II. NONSMOOTH DYNAMICS - THE CONVEX 
APPROACH 


Consider a piecewise smooth dynamical system 

^ _ / f+i^) for > 0 ’ ni 

dt ^ \ f-{x) for v{x) < 0 . ^ 

This paper is about how we extend such a model across 
v{x) = 0, and the issues that arise in doing so. 

In the standard Filippov approach prominent in 
variable structure control [1, 12^ 1^ and in piecewise 
smooth dynamical systems theory ^ Q , when a system 
switches between two systems as in o, we form their 
convex combination 


dx ,, f+{x)+f_{x) , f+{x)-f_{x)^ 

It = rfeA) = -^^-A 

( 2 ) 

(or a similar form ^ = [/+ (x) + /_ (a:)] -f 

u [f+ix) — f-ix)] with X = 2u — 1), where 


J sign (v) it V 0 , 
{ [-I,+I] if n = 0. 


( 3 ) 


thus 


/(a;; -kl) = f+{x) and f{x;-l) = f-{x) . (4) 


The standard approach then seeks so-called sliding modes 
which satisfy = 0 on n = 0. Taking coordinates 
X = (xi, X2, ■■■, Xn) such that a:i = v, and writing / = 
(/i, / 2 ,..., /rt), the dynamics of sliding modes is therefore 
given by 


0 = h{x-X^) 

dj{x 2 ,.:,Xn) = (/ 2 (a:; A®), ...,/„(a;; A®)) 

( 5 ) 

for some A'* S [—1,4-1]. If no solutions of the sliding 
mode problem exist for A® S [-1,4-1], then the flow of 
@ crosses through v = 0 transversally. 

This definition of sliding and crossing at a switching 
surface is sometimes justified by smoothing (or ‘regular¬ 
izing’) the combination by replacing A with a smooth 
sigmoid function, and showing that slow dynamics of the 
smoothed system contains invariant manifold dynamics 
equivalent to the sliding mode dynamics ([S]), (for theo¬ 
rems regarding this equivalence, which remains impor¬ 
tant in spite of what follows, see [H, [2^, H^). Unfor¬ 
tunately such an argument is tautologous. Smoothing 
simply preserves the dynamics we have imposed already 


on xi = 0 


by writing the expression (ED, and we shall see that by 
changing Q we can not only obtain different dynamics, 
but that the difference persists under smoothing. Two 
examples illustrate the ambiguities of smoothing a dis¬ 
continuity very clearly, and while almost trivially simple, 
they help motivate the general arguments that follow. 

Example 1. Let us start with a smooth system that we 
wish to study, 

j-^{xi,X2) = {ipeixi),l-2if‘l{xi)) , ( 6 ) 

where e is a small positive parameter and ipg(xi) is a 
sigmoid function sueh that (po^xi) = sign(a;i) (so thee —>■ 
0 limit of ipg[xi) is ^). This has an invariant line xi = 
ec, where c is a constant such that ipg{sc) = 0, on which 
the dynamics is given by 

^(xi,X2) = (0,1) . (7) 

This is illustrated in the far right of figure]^ Now con¬ 
sider a piecewise smooth model of In the limit e = 0 
equation becomes 

|j(a;i,X 2 ) = (sign(a;i),-I) . (8) 

Filippov’s linear combination for (0), using simplifies 
to 


^(xi,X2) = (A,-l) , (9) 

where A is given by 0). A sliding mode on the discon¬ 
tinuity surface Xi = 0 must satisfy ^xi = 0 by 0), so 
sliding modes exist for this system with A = 0, giving a 
sliding vector field ^{xi,X 2 ) = (0, —1), shown in the far 
left of figure [H Clearly this is not even close to being 
similar to the smooth system’s dynamics in the far right 
portrait in figurel^ given by m- Moreover if we smooth 
this system, by replacing A with a sigmoid function (we 
can use (fg again for brevity), we obtain 

^(xi,X2) = ((Pe(xi),-1) , (10) 

a system with an invariant line at some xi = c, where c 
is a constant such that ips{xic) = 0, but on this line the 
vector field is again ^{xi,X 2 ) = (0,-1), inconsistent 
with &; this is the second phase portrait in figure[J\ 

We obtained inconsistent sliding dynamics here, which 
persists under smoothing, because we did not keep the 
hidden term ip')., neglecting is because is just unity for all 
xi If we take (0) again, and carefully replace (pe{xi) 

with A in the limit e —>■ 0, we obtain 

^(xi,X2) = (A,I-2A") . (11) 

A sliding mode on the switching surface xi = 0 must 
satisfy ^xi = 0, and this has solutions for A = 0, which 
now gives a sliding vector field jg{xi,X 2 ) = (0,-1). This 
correctly captures the dynamics of the smooth system 0 
at the transition, and is illustrated by the third portrait 
in figure [7J Moreover, obviously if we now smooth this 
by replacing A with a smooth sigmoid function, we regain 
(0), the far right portrait in figure [H 
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nonlinear: 


from which the presence of a discontinuity is not even 
evident. Filippov theory therefore does not apply, and 
the system crosses through xi = 0. If we more carefully 
replace (pe by X in m, we obtain 

^(a:i,X2) = (2A2-l,l), (14) 

which has two sliding modes (where ^xx = Oj given by 
X = implying ^{xi,X 2 ) = (0,1), equivalent to 

the invariant dynamics of the smooth system. Moreover, 
when this system is smoothed by replacing X with a sig¬ 
moid function like ife, we of course return to the original 
system m- 



FIG. 1. The system ^(xi,X 2 ) = (sign(a:i), 1 — 2(sign(xi))^) 
viewed ignoring the nonlinear dependence on sign(xi) (left) or 
respecting it (right), in each case showing the discontinuous sys¬ 
tem and its smoothing. The phase portraits as shown correspond 
to equations ([9|, ()lQp . ()lip . ([6|, respectively. 


In this example we obtained sliding dynamics that was 
not only quantitatively wrong compared to the system we 
were trying to model, but traveled in entirely the wrong 
direction. This happened because we neglected terms 
nonlinear in the switching term or X. 

The point is that if we start from the piecewise smooth 
model (IH), then either ® or (El) are valid ways of con¬ 
tinuing the system across xi = 0, amongst an infinity of 
other choices. In the following sections we will show how 
to express those possible systems, separate ‘linear’ and 
‘nonlinear’ effects, and we shall highlight some more and 
less subtle effects they give rise to. 

A similar example to that above reveals that nonlin¬ 
ear terms can even affect whether or not the switching 
surface is crossed. 

Example 2. Consider a system 


^(X1,X2) = (2v5^(xi) - 1,1) . (12) 

The idea is similar to the previous example so we just 
outline the calculations. This system has two invariant 
sets xi = sc± such that </5e(ec±) = ±l/-\/2, on which the 
vector field is ^{xi,X 2 ) = (0,1). In the limit e = 0 we 
have 


(13) 


linear: 


Eg2. 


smooth 












0 

-E 6 +£ 


nonlinear: 



-T- >X 

0 



FIG. 2. The system ^(xi,X 2 ) = (2(sign(xi))'^ — 1,1) viewed 
ignoring the nonlinear dependence on sign(xi) (left) or respect¬ 
ing it (right), in each case showing the discontinuous system 
and its smoothing. The phase portraits as shown correspond to 
equations II13D . II13I) . (1141) . (1121) . respectively. 


In this example the piecewise smooth model outside 
xi = 0 is given by which is formed by patching 

together two systems for Xi > 0 and xi < 0. Even 
though the vector fields appear continuous, their regions 
of definition are not. We may then assume the system 
is continuous, or consider dill) (and again an infinity of 
other choices) as a valid continuation of the system across 
xi — 0. In this case, the assumption of simple continuity 
turns out to be incorrect for modeling the system ED 
that we started from. 

These are just trivial systems, and in the next section 
we will show how vast the field of possible regulariza¬ 
tions is. Yet the ambiguities these examples highlight 


^(xi,a;2) = (1,1) 
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are extremely important to understanding the modeling 
of switches. Other examples can be much more subtle, 
as we shall see later. In section IIVI we discuss how to 
handle nonlinearities to avoid ambiguity. 


III. THE SMOOTHING TAUTOLOGY - 
KEEPING THE HIDDEN TERMS 


What is not clear in many works concerning smoothing 
is that there are many ways a system can be smoothed, 
and that they result in non-equivalent systems, any of 
which is equally valid as an approximation of the piece- 
wise model ©■ In this section we make this ambiguity 
explicit in the piecewise smooth system, as a ‘hidden’ 
term in the vector field. 

Consider smooth systems expressible in the form 


^ = f s.t. 




sign (u) if |u| > e 

[-1, -fl] if juj < £ 


-I- 0 {e/v,e) 


(15) 


where the transition function is at least continuous, 
and e > 0 is small. The error term 0{e/v,e) ensures 
that we obtain o as £ —>■ 0. We do not specify the form 
of (pe beyond this, for example it need not be monotonic, 
and need not be bounded to [—1, -1-1] at u = 0. 

Let us moreover suppose / to depend on in an arbi¬ 
trary manner, that is, it could be a polynomial, trigono¬ 
metric, exponential or other function of the sigmoid. 
Nevertheless, if we take the limit £ —?> 0 in m we ob¬ 
tain ©■ Does the dynamics of (fT51) therefore tend to the 
dynamics of © as described by @? 

Since we do not know precisely the form of in the 
limit £ —> 0, let us introduce a simpler ‘reference’ sigmoid 
for which we do. We define it as 


Ae{v) 

S.t. 


J sign(t?) if |i;| > e 
\ [-1,+1] if \v\ < £ 

A'^{v) > 0 for |u| < £ , 


+ 0 (^) j 


(16) 


and we shall use this to analyze the behaviour of pe as e 
tends to zero. Since A£(u) is differentiable and monotonic 
for £ > 0, at least for |u| < £, it has an inverse I^(A£) 
such that V{Ai;{v)) = v/e. We can then define a new 
expression for as a function of Ag, 


4-(Ag) := (£U(Ag)) . (17) 


Let us now use the reference sigmoid Ag to form a linear 
combination of f{x; ±1), 


/:(:r,Ag(u)) = i±4iM/(^;+i) + L4£M/(a,;-l)+... , 


taking the difference between this system and dig , we 
have 


^ ^ C{x,A,{v)) + £ {x,A,{v{x))) , (18) 

where £ {x, Ag) = / (a;; 'L(Ag)) — £(x, ±1) vanishes when 

Ag = ±1. 

Because the reference sigmoid Ag varies monotonically 
with V and is well-behaved as £ tends to zero, we can 
identify it with a new variable A, then dT51) becomes 

^ -M) -f ^-^/(x; -1) £ (x. A) . (19) 

The quantity A will vary dynamically according to ^ A = 
where ^ is strictly positive (which would not 
have been the case if we had identified A with p^). 

The expression (1191) is now e independent, and we may 
consider A to represent any p^, including the limiting 
function A = (j)o{v) = sign(u), with which (fTOl) defines a 
piecewise smooth system consistent with ©, but notably 
more general than ([2]). 

Since £ vanishes for u 0, it is not fixed by (P). f is 
the hidden part of the vector field, and by changing it we 
obtain arbitrarily many different non-equivalent dynam¬ 
ical systems, all consistent with ©. 

There is of course no contradiction or paradox, nor in¬ 
deed any surprise, in the fact that infinitely many differ¬ 
ent smooth systems expressible as m, have the same 
limit © as £ —>■ 0. This is merely due to the non¬ 
uniqueness of the singular limits by which smooth func¬ 
tions may tend towards discontinuities. What is impor¬ 
tant is that we can represent these different systems in the 
piecewise smooth limit by dig. The system (P, there¬ 
fore, is only the simplest member of a general class given 
by (dg. 

The argument above began with (1151) as a prototype 
smooth system, but one may consider many other kinds. 
For instance consider a number of switches that activate 
at a common threshold x = 0, with each switch having 
a different sigmoid-like behaviour pi,p 2 ,..., so (fT51) be¬ 
comes 

fl = f ix;piivix)) ,p 2 {v{x)), where 


‘Psk (x) e 


sign(x) if |x| > £fc 

if \v\ < £k 


For example we might let 


+ 0{ek/v) . 

( 20 ) 


Psi (v) 

(x) 

Psaix) 

‘Pe4(v) 


1 , 

2 ' 2 l-|-(j)-l)l/ei ’ 

tanh(x/£2) , 

Erf(x/£3) , 

§ arctan(v/£ 4 ) , etc. 


We can again use a reference sigmoid Ag with inverse 
V = eV{Ae) to define a new expression for each t/jg^. 


such that ^ = £(x,±l) in the limit |x|/£ —>■ oo. This 
will not be exactly equal to m for general v and £, so 


4'fe(Ag) := p^^ (£E(Ag)) 


( 21 ) 
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and derive dun as above with £{x,A^) = 

/{x;^l(Ae),«'2(Ae),...) - (l+Ae)/(x; + l)-(l-A.)/(.;-l) ^ 

which vanishes when Ag = ±1, resulting again in the 
piecewise smooth system (HU). 

An interesting example arises when the sigmoids ap¬ 
pear only linearly, but in a matrix form such as 

A = +1, + _1) . 

( 22 ) 

which we may call a matrix regularization of o, where 

M{v) = 

or even 

M{v) = 


[ {v) 

0 

0 


0 

<7>e2 (v) 

0 


0 

0 



V : 




/ (^) 

<5£12(^) 

Xisiv) 


^£21 (^^) 

(v) 



deal (^') 

X 32 i^’) 

Veu (v) 


V : 





where the off-diagonal elements are gaussian-like func¬ 
tions Sij{v) = 0(e) for I'cj > e and 6ij{v) £ [0,1] for 
I'cj < e. To handle these off-diagonal terms we need an 
extra step, because they are not sigmoids, nevertheless 
we can recast them again in terms of the reference sig¬ 
moid Ag via a function Ay := dy {eijV{Ag^)), then to 
derive the corresponding hidden terms in £" is a simple 
exercise, resulting in (fTO|) as before. 

The class of piecewise smooth dynamical systems that 
arise from such considerations is clearly vast, but all cases 
are encompassed in (HU with a suitable hidden term 
£ [^ . Of course above we have written £ in terms of 
the systems (HSl), o, or (1^ . expressions that are not 
useful in themselves (except that they show £ vanishes 
outside the switching surface). Some means to express 
the hidden term £ explicitly is now required and, where 
possible, to determine it from a system’s behaviour. We 
broach this in the following section. 


IV. SIGMOID ASYMPTOTICS 

Given a system ^x = f, and from only the knowledge 
that / jumps between functional forms f+{x) and f-{x) 
as v{x) changes sign, we wish to derive a complete model 
for the system. This will take the form / = /(cc; A) in 
terms of a piecewise constant A, or / = f{x;(pg{v)) in 
terms of a sigmoid function ifg. 

Below we will attempt formally to expand / as a power 
series 

, OO 

-£-= f{x-,X) = '^an{x)X^{v{x)) (23) 


where for a piecewise smooth model we assume A = 
sign(w), and for a smooth sigmoid model we replace 
A I—>■ (feiv), assuming ipg{v) > 0 and 

</3e(z;) = sign(t!)-I-0 (e) for |t!| > e . (24) 

The possibility of approximating a function of several 
variables {x, y) by a sum of functions of a single variable 
v{x) has arisen particularly in the context of universal 
approximation by neural networks [3, B ) where the 
sums of interest include sigmoid functions ip{vi{x)) that 
may each have a different threshold Vi = ki ■ {x — Xi) 
and stiffness [fc^j for some vectors 2 (falling 

therefore into our subclass of matrix regularizations in 
(l22)l b Here we are interested in a single switching thresh¬ 
old, say Xi = 0, and will be interested in the limit of in¬ 
finite stiffness \ki\ —^ 00 . An alternative is to use a series 
of powers of the sigmoid function (/j(u), and the use of 
such series as solutions of nonlinear differential equations 
has been explored in [2l|, for ip = tanh. 

The task now is to find the coefficients of the power 
series (l23ll . 


A. A sigmoid sum 

Comparing (|Mll to m in the states A = ±1, we have 
immediately that 

00 

a-nix) = f+{x) (from v > 0) , 

(25) 

a„(a;)(-l)" = f-{x) (from z; < 0) . 

n—0 

The sum and difference of these give 2 J2n even C(n = f+ + 
/_ and 2 odd On = /+ — /_, allowing us to eliminate 
ao and oi to obtain 

^ = i(/+ + /_) + i(/+ - /_)A(n) + £ix; A) , (26) 
where 


£{x-, A) = (a 2 n + X{v)a 2 n+i) (</5^"'(w) - 1 ) 

n—1 

00 n — 1 

= (a^-i)EE [a 2 n{x) + Aa 2 „+i(a;)] 

n—1j=0 

:={X^-l)g{x,X) , (27) 

simply factoring out A^ — 1 in the second line, and defin¬ 
ing the double sum as a function g. Note that while £ 
vanishes for A = ±1, i.e. outside the switching surface, 
the function g can take any value. Hence we obtain ex¬ 
actly (HU), and a formula for £. 

In the piecewise smooth case, the remaining coeffi¬ 
cients a „>2 that make up the hidden term g may be 
determined if we have any information about f at v = Q. 
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For example, say it is known that for A = 0 we have 
/(2;;0)|„^o = r{x), then 03 = i(/+ + /_) - r(x), and 
successive coefficients can be eliminated with knowledge 
of derivatives of / with respect to A. If these quantities 
cannot be derermined by direct observations of a system, 
we may propose forms of g to fit dynamical observations 
or physical laws, depending on the modeling problem. In 
the smooth case, where A is replaced by a differentiable 
sigmoid function then more can be found by consid¬ 
ering the asymptotics of the transition. 


B. Asymptotics of the series 


Here we consider how (ESD-dUl) may be applied to a 
system where the transition between f± is smooth. 

In a system that makes a sudden transition between 
behaviours ^ = f+{x) for large positive g = v/e, and 
^ = f-{x) for large negative g = v/e, let us assume 
that /+ and /_ are actually steady states f{x;y+) and 
f{x] u-) of a larger bistable system in x and y, such that 


y 


?/+ as 77—^ -koo , 
y- as ?7 —?> —00 , 




dx 

dt 


/+ as ry —>■ -koo , 
/_ as ry —?> —00 . 


( 28 ) 


The vector y may represent microscopic properties or ex¬ 
ternal forces independent of x, which vary significantly 
only where their value jumps. The scalar g may merely 
be a fast timescale t/e for small e, or else some function 
g(x, e) of the internal state. Without specific knowledge 
of the y system, we seek to model the effect on the system 

t = m- 

If we could linearize with respect to y near the steady 
state /+ = f{x]y+) for large negative g, treating g 
as a time-like variable, then / should diverge from /+ 
along a direction in phase space with a time con¬ 
stant K (the Lyapunov exponent). The resulting approx¬ 
imation df /dg ^ — f+) has a solution of the form 

/ -^ /+ -k 9+e'^'^ l^rig~'^ for some coefficients /3„. 
This is easily extended if we cannot linearize. For exam¬ 
ple, if / = /+ is a degenerate steady state of the y system, 
the same asymptotic formula applies with k = 0 if / = /+ 
is a saddle-node, where df /dg ^ \ f — f+\{f — f+)/ldi for 
some constant /3i. The same formula again holds with 
the exponent replaced by —K\g\P near a non-hyperbolic 
steady state of the form df /dg ^ K\g\P~^{f — /+). 

Based on the typical divergence from a linear or non¬ 
linear steady state, therefore, near either of the states 
^ = /+ or ^ = /_, the system behaves asymptotically 
like 


dx 

dt 


Mx)+g^ix)e ^(^\x){e/v) 

n—0 


( 29 ) 


as u/e —>■ —00 for i = -k or u/e —>■ -koo for i = —, for pos¬ 
itive constants Ki,pi,[ 3 n'^. The vector-valued functions 


gi{x) describe the trajectories in function space by which 
^ departs from each of the fi. Provided /3q*^ = 0 if 
Ki = 0, this ensures that the second term vanishes for 
v/e —>■ ±00. 

If we had the exact differential equations for f{x',y), 
then we would now apply asymptotic matching (see e.g. 
Q) to unify the two approximations in (j29ll . When no 
such differential equations are known, we may seek to 
represent / by a series approximation in x only, since y is 
assumed to be beyond ready observability. The terms of 
the series must be functions that give the correct asymp¬ 
totic behaviour for positive and negative v, and this can 
be achieved with sigmoid functions tp that vary monoton- 
ically from —1 to -kl, switching value at g{x) = 0. Thus 
we take the series (j23p with (l24)l . 

Finding the first two coefficients is similar to the piece- 
wise smooth case. Comparing (l23t to (l29l) at the steady 
states /+ and /_ we have that oq = ^{f+ + /-) + X) 
and ai = i(/+ - /-) + X) a 2 n+i, so 

^ + ^^—^^<fe{v)+£{x,(pe{v)) (30) 

where S{x,(pe) = {(fi^ — l)g{x,(pe), with E and g being 
the same functions as in the piecewise smooth case (EH). 

The remaining coefficients gathered together in g can 
be found by assuming that has an asymptotic expan¬ 
sion of the form 


ipeix) ^ sign(z;) ( 1 -k e 7n(e/'c) 


n—Q 


( 31 ) 


consistent with (IMl) . For example, the Hill functions 
[III behave as Z{x) = ^ife+gi/e ^ 5 + sign(a; - 
6) (i - -k ...) for a; > 0, and we let v = x-6, 

fe(v) = 2Z(0e^^‘^) — 1 , and other common sigmoids in¬ 
clude (ps{v) = tanh(u/£) ^ sign(u) (l — -k ...), 

<Pe(u) = Erf(w/£) ^ sign(w) — ((e/u) — ^{e/v)^ + 

...)/V^. 

To give an example of the coefficients an> 2 , if we as¬ 
sume some symmetry in the asymptotic limit such that 
K_i_ = K_ and y 7 + = p-, then it makes sense to choose 
Kj = K+ = K- := K and p^ = p+ = := p, then com¬ 
paring terms of order allows us to find 02 and 03 

as 

a2 = 4^ + 9 -Po^') - Em=2 , 

«3 = 4 ( 5 +/?^“^ - ) (32) 

-\{f+ - f-) - YZ=2^0i2m+l , 


and so on, recursively for each coefficient 


V. EXAMPLE: SWITCHING ASYMPTOTICS IN 
A CONTROL CIRCUIT 

The most obvious kind of sigmoidal behaviour in 
physics is a mechanical or electronic binary relay. Con- 
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sider a typical electrical model based on a DC-DC con¬ 
verter, with a battery voltage Vf, and inductor current 
I. A connecting circuit, with a capacitance C and resis¬ 
tance R across which the voltage is V, is disconnected if 
V exceeds a value 14. The empirical model of the on and 
off positions is 




on 


“off”: 


r L^j = Vo-V, 
\RCiv = IR-V, 


/ Lf-J = Vo , 

I Rciv = -V■ 


(33) 


becomes (l35|) as £ —>■ 0, i.e. this is an ^-perturbation of 

(ESI- 

The system (IM)) has a small deviation from the steady 
states n — 0 and ^ = 1 for |u| > e, described by the 
derivative with respect to V/e, 

^/=(-V',/Wrt)5^ + 0(.), (37) 

The p' term means this deviation will be a dependent. 
The term means the deviation will be of order 

0(|14 —td|/e) or smaller, as discussed in the previous 
section leading to (l2^ . Take for example (ps to be an 
arctan function, specifically 



ifeiVb -y) = k + ^ arctan 

= step(14 - 1^) - + 0 

then 

dp _ dipejVb - V) _ 7r~^ _ 

dV/e ~ dV/e ~ 1 + \Vb - V\^l£^ 




0 


{(v,-vr) > 


What then happens during the transition between 
these two? The simplest model we might write down, 
and the canoncial Utkin/Filippov approach [1, ex¬ 
presses (1^ as 

{LjjI,RCj^V)=f = iVo-V) + {-V,IR)p 

(34) 

where p = step(14 — V) ■ 

The switch p is related to the multiplier A in the previous 
section hy p = ^ + ^X, we use p here as it is more common 
in this application. 

Let us then consider a more general class of systems 

{Lj-J,RCj-^V) =f = {Vo,-V) + i-VpJRpip)) 
where p = step (14 — V) , 

(35) 

such that p{p) = step (14 — V)- (This class is introduced 
for demonstration, and it may be possible to motivate 
more experimentally relevant classes from the precise 
physics of a given device). For the sake of an explicit 
example we shall consider p{p) = p — (t(1 — p)p below, 
so that the original system (IMl) corresponds to ct = 0 . 
We shall ask to what extent we can distinguish between 
systems with different a. 

The smooth function replacing the step should be 
monotonic and continuous, say p = (psiVb — V) defined in 
terms of a stiffness parameter e, such that v?o(14 — Id) = 
step (14 — y)- The resulting system 

{lI-j,rc!Lv) = f = ( 14 ,- 1 ") + {-ypjRp{p)) 

where p = ipe{yb-y), 

(36) 


hence 

QYjJ ~ (^)) 2\Vb-vp T® ((V6-y)3>^) ■ 

(38) 

The difference between models with different a is there- 

2 

fore very small, shrinking with e and with for ^ 

away from 14 (and alternative sigmoid functions may give 
deviations exponentially rather than polynomially small 
in |14 — l"|/e). Observing the deviation directly, there¬ 
fore, would be a challenge in practice. Nevertheless we 
will show that a has significant and observable effects on 
the stability of the circuit no matter how small e is (i.e. 

how close the system is to discontinuous), because the 
2 

term permits significant deviations near V = Vb. 

So we now turn to the circuit’s behaviour. 

With the switch in the “on” position the state spirals 
in towards an attracting focus at / = Vq/R, V = 14, as 
shown in figure 01 In the “off” position the system is 
driven towards the switching surface 1" = 14. We then 
need a law describing how the systems evolves at the 
switch. 

We can use the smooth system to derive how p jumps 
between 0 and 1 across the switching surface. Differenti¬ 
ating p = (ps{Vb — V) with respect to t gives 

d dV d d j 

where e = satisfies e < e < 4e inside the transi¬ 

tion region |ld — 14| < £, hence e = 0 (e) is small around 
the switch. We can therefore introduce a fast timescale 
r = t/e in terms of which ^ = ^ = / • VI". The de¬ 
pendence on s is supressed (but for e = 0 the r timescale 
































is instantaneous compared to t), and we can take the 
limit e —>■ 0 to obtain the piecewise smooth system, in 
which the transition region becomes the switching sur¬ 
face V = 14, on which we have 

= {Vo,-Vt) + {-Vb^l,IRp{^l)) (40) 

Thus (HDl) prescribes the dynamics on (/, p) at a point 
on the switching surface V = Vb. Let us take the function 
p{p.) as 

M (41) 

as an example, and assume |(t| < 1. Then for IR <Vb we 
have dp/dr < 0, implying that the flow crosses directly 
from V > Vh to y <14- For IR > Vt instead, the 
interval or layer of the transition, p € [0jl]i contains 
an attractor in the form of a saddlepoint at /r = ^, 

IR = VoVt-crVoiVt+Vo) - dynamics outside the layer 
are shown in the main part of figure lU with a blow up 
(inset) of the dynamics on p G [0,1] inside the switching 
surface V = Vb, shown near the saddle. 

Solutions that enter the layer p G [0,1] will be at¬ 
tracted to one of the two unstable manifolds of the sad¬ 
dle, either traveling toward increasing / with V remain¬ 
ing on the switching threshold (shown in figure SDa)), or 
traveling toward decreasing I, before detaching from the 
switching surface at IR = Vb and spiralling in towards 
the focus in the “on” position (shown in figure |4l[b)). 

The critical effect of a in this system is to shift the cur- 

rent I = y^y^_J'y^^y,,+Vo) the saddlepoint occurs 

on the switching surface. Figure U shows the behaviour 
of the system from an initial condition {I,V) = (0,0) for 
two different values of a. 

If we set cr = 0, then from I = V = 0 the system tra¬ 
jectory enters the sliding mode to the right of the saddle, 
and the current grows unboundedly, while the voltage re¬ 
mains fixed at y = 14 . In this case the switch depends 
linearly on p, and the sliding mode is simply as would be 
obtained from Filippov/Utkin’s standard method [§.1^. 

where the saddlepoint at p = -^,IR=^ is usually 
called a psuedoequilibrium. 

For cr nonzero, instead, the saddlepoint moves to a 
higher current, and from the initial condition I = V = 0 
the system hits the sliding mode to the left of the saddle, 
the current falls until the sliding mode terminates, and 
the system then spirals in towards the focus. The differ¬ 
ence in the dynamics at y = yj, has to do with the fact 
that, with a Q, the system is nonlinear in the switching 
parameter p. 

The effect of nonlinearity in the voltage response to the 
switch here is relatively minor - shifting the saddlepoint 
- but with stark effects on the system stability. The ef¬ 
fects can easily be more interesting, as we show with an 
abstract but novel example in section IVlIl First we shall 
generalize the methods above for analyzing the piecewise 
smooth system, without having to concern ourself with 




FIG. 4. Phase portrait of the current, showing the piecewise 
smooth system (main picture), and a blow up of the switching 
surface into the layer p G [0, 1] revealing a saddle. The trajectory 
from / = y = 0 is shown, simulating values L = 5, Vb = 5, 
Vb = 6, RC = 5/2, R — 15/4, and (a) <7 = 0, (b) a = 1/2. 

the asymptotics of smoothing, yet without losing the dis¬ 
tinction between systems that behave differently near the 
switch. 


VI. NONSMOOTH DYNAMICS - THE 
NONLINEAR APPROACH 

We can generalize the methods used in the previous 
section with a slight development of methods proposed 

in [1, [T^. 

We have shown in section m that we can express the 
system o as 

dt 2 2 \ ^ J \ J 

where 

A = sign?; and v£{x] A) = 0 . (43) 

The term £ is the ‘hidden’ part of the vector field, the 
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condition v£ = 0 meaning that this vanishes outside the 
switching surface. 

At u = 0 the system o with ([3|) gives a differential in¬ 
clusion, but (as we found in the last section) it turns out 
that the value of A can be fixed in most cases by dynam¬ 
ical arguments. For £ = 0 this is done using Filippov’s 
sliding modes and in these ideas were extended 
to the case when £ is nonzero, i.e. when depends 
nonlinearly on A. We shall review the method briefly. 

We introduce a transition timescale t which is instan¬ 
taneous in the original timescale t (so we may think of 
this as r = lim£_>.o t/e for some e > 0 , as we obtained in 
the previous section). Since A switches between ±1 as u 
changes sign, its variation is driven by the component of 
the vector field along V?; normal to the switching surface, 
so we define the dynamics of A as 

-^A = f(x; A) • Vv on ?; = 0 . (44) 

dr 

This is the generalization of the expression we obtained 
in the circuit model in section |Vl and can be more gener¬ 
ally derived from a smoothing of © using arguments of 
asymptotics and geometric singular perturbation theory; 

see [I 1 -II 3 - 

The expression (1421) with (l44l) is now sufficient to 
prescibe the dynamics of the piecewise smooth system 
fully. Combining them at the switching surface, taking 
coordinates x = {xi,X 2 , ■■■,Xn) in which xi = v{x), and 
letting / = (/i, / 2 , ...,/„), we have the two timescale sys¬ 
tem 

j-^{x 2 ,:;Xn) = (/ 2 (a:;A),...,/„(a;;A)) J ^ 

on Xi = 0. When a solution reaches the discontinuity, 
the system (l45l) facilitates either the transition from one 
side of the switch to the other, or, if equilibria of the r 
timescale subsystem exist (i.e. where ^A = 0), then A 
collapses to a value A'* given by 

0 = h{x;Xn \ , . 

^(x2,...,x„) = (/2(a:; A"), ...,/„(a;; A®)) J 

on xi = 0. These equilibria, if they exist, represent slid¬ 
ing modes, whereby the dynamics sticks to the switching 
threshold v = 0 and evolves {slides) along it. We see 
that these are precisely the sliding modes defined in ([5]) 
by Filippov, except that now with dH]) and (H51) they 
apply also to systems with nonlinear dependence on A. 

These are the basic elements for solving a system with 
nonlinear switching. We have seen in three examples how 
nonlinear terms can alter the dynamics. Nonlinearity can 
head to much more interesting dynamical phenomena, an 
example of which is given below. 


VII. EXAMPLE: HIDDEN DUEFING 
OSCILLATOR 

Consider the system 

^{xi,X 2 ) = {x 2 - cxi, -X^ - bx 2 + acost) , (47) 

where A = sign(xi), for constants a,c,b. The dynamics 
is that of a fused focus, that is, a folding of the dynamics 
either side of the switch that creates focus-like attraction 
towards Xi = X 2 = 0 ; this is illustrated in figure [5Ki). 

For xi ^ 0 this is indistinguishable from the linear 
switching system 

^{xi,X 2 ) = ix 2 - cxi, -X- bx 2 + acost) , (48) 

where A = sign(a:i), which is well described by Filip¬ 
pov/Utkin sliding mode theory [iM 1^. 

Both (l47l) and (l48l) are clearly models consistent with 
the system 

(j / '1 — / (^2 — cxi, —1 — bx 2 + acost) if xi > 0 , 

dt^xi,X 2 ) y {x 2 — cxi, -\-l — bx 2 + acost) if CCi < 0 . 

(49) 

Let us apply the methods from the previous section to 
(iTfl) to see the effect of the A^ term in (iTfl) . Substituting 
(SZl) into (HSI) to find the dynamics on xi =0 gives 

{■j-X, ^X 2 ) = {x 2 , —X^ — bx 2 + acost) on xi = 0 , 

(50) 

where r is an instantaneous transition timescale. Only 
the attracting point xi = X 2 = 0 satisfies the equa¬ 
tion (H51) , so sliding modes are confined to the origin and 
would seem to be of little interest. However, in the vari¬ 
ables (A, CC 2 ) the system (l50l) is a form of Duffing oscillator 
(see e.g. [l 2 |), and this can lead to rich dynamics inside 
the transition layer A G [—1, -fl]. 

Simulations indeed reveal nontrivial dynamics on A 
and X 2 for A S [—1,4-1]. One simulation is shown in 
figure EKii), showing the switching surface xi =Q ‘blown 
up’ into a switching layer A G [—1,4-1] on xi = 0. The 
variable X 2 varies only slightly, while A oscillates in an 
irregular pattern between values of around JAj < a^/^. 

The structure of the orbit in figure [5Kii) is not visible, 
and may appear to be of little significance anyway, since 
it is negligible in the state space of {xi,X 2 ) in figure [SKi). 
The switching paramater A may, however, have a notica- 
ble effect on the system, for example if it affects another 
variable. Let us introduce a third variable satisfying some 
simple equation involving the switch, say 

= X- X 3 , (51) 

where /r is some constant, then 0:3 will be simply 4il 
outside the switching surface, but inside the surface it 
will track the dynamics of A. We choose /i small enough 
in our simulations so that X 3 and A are indistinguishable. 

After discarding transients, 2:3 (t) is seen to have oscil¬ 
latory, complex, and chaotic dynamics for different pa¬ 
rameters, including Ueda’s chaotic attractor [ 2 ^, as ex¬ 
pected from a Duffing oscillator. FigurelHKa) shows a plot 
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FIG. 5. The flow of: (i) the fused focus (|47p . and (ii) its switch¬ 
ing layer II50I) . shown at a fixed time, for constants c = 0.1, 
b = 0.05, a = 0.15. An orbit is simulated over a time t from 0 
to 500. 


of ^3 (f) . The speed of the timescale t matters here, so we 
let T = t/e, the ideal instantaneous switch corresponding 
to the limit e —f 0. The curves (i) and (ii) show X 3 {t) (or 
equivalently the value of A) for e = 10“^ and e = 10“®, 
respectively. The oscillations become more regular as e 
tends to zero, with the amplitude of this regular oscilla¬ 
tion tending to around a}^^. 

Contrast this with the curve (iii), which shows the 
graph of xz{t) for the linear model (H51) . obtained by re¬ 
placing with A in O, and gives instead always a 
steady oscillation (which does not vary appreciably with 
e), but with an amplude a, rather than the of the 
nonlinear system. 

Throughout this paper we have tried to show that 
smoothing is a tautologous process, that is, it cannot be 
used to justify one way or another of resolving a discon¬ 
tinuous system, say taking (SZl) or (H51) to resolve (ggi) . 
and that systems must be distinguished in the expres¬ 
sion of the discontinuous systems themselves. We demon- 




FIG. 6. Simulations of hidden dynamics, plotting z{t) (or 
ip{x{t)/e) for: (i) the nonlinear system with e = 10~^, (ii) 
the nonlinear system with e = 10®, and (i) the linear system 
the e = 10“^. Panel (a) is a simulation of the switching layer 
system (I50II . and panel (b) is a simulation of the full system 
(1471) smoothed out using A = ipe{xi). All simulations are with 
c = 0.1, h = 0.05, a = 0.15. 


strate this finally by smoothing out the systems (HZD and 
(H51) directly, replacing A in those expressions with 



sign(a;i) if |a;i|>e, 
xi/e if \xi\ < e , 


(52) 


and we simulate in figure Efb) the orbits correspond¬ 
ing to those in figure Ha). The qualitative dynamics 
is similar, with only minor quantitative differences in 
the precise form of smaller oscillations, confirming that 
smoothing only preserves the differences between the lin¬ 
ear and nonlinear models. Simulations using alternative 
smoothing functions such as (ps{xi) — ^ arctan(xi/£) or 
(pe{xi) = tanh(xi/£) yield results with no significant dif¬ 
ference to those using (15^ . 


VIII. CLOSING REMARKS 

The traditional use of Filippov’s theories to tackle 
piecewise smooth systems neglects the possibility of non¬ 
linear dependence on A; (we should remark that Filip¬ 
pov’s work 0 sets out results for much more general dif¬ 
ferential inclusions, but most of his detailed theory and 
most modern theory derived from it assume linear depen¬ 
dence). In many cases the assumption of linear depen- 
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dence on A will be sufficient, but nonlinear dependence 
may have to be introduced if a linear model’s dynamics 
that is inconsistent with experiment, if we smooth the 
system out, or, something we have not discussed here, 
if singularities arise that exhibit degeneracy inside the 
switching surface [l^ . 

Great intricacy is to be found in not only the global, 
but the local, effects of discontinuity, due to singularities 
[T^ , bifurcations [1, [13, Ull , and now it seems also due 
to the hidden dynamics of the switch itself. To analyse 
these models and apply them to physical systems requires 
looking closer at how piecewise smooth models are formu¬ 
lated. The regularizations in section Hill indicate that the 
range of models possible is much wider than sometimes 
appreciated, while the sigmoid series expansions in sec¬ 
tion IIVI attempt to form a general approach to encoding 
these different systems via nonlinearity. 

In conclusion, whether we model a system as smooth 
or nonsmooth is not the crucial question, as was illus¬ 


trated in section [Vn] in particular. What matters is that 
hidden terms - those which are negligible away from the 
switching surface - are respected in whichever modeling 
framework we choose. Section llVI provides the means to 
preserve these hidden terms in a piecewise smooth model, 
and to discover them from small deviations in the ‘real’ 
system if it is actually smooth. If we choose to analyze 
a system using piecewise smooth theory, then section lyTl 
outlines the basic approach, and the examples through¬ 
out this paper are merely a glimpse of the kind of strange 
behaviours that remain to be explored, and which may 
perhaps offer new means to model the dynamics of irreg¬ 
ular systems. 
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